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Abstract 

We investigate numerically and analytically Potts models on "thin" random graphs - generic 
Feynman diagrams, using the idea that such models may be expressed as the N — * 1 limit of a matrix 
model. The thin random graphs in this limit are locally tree-like, in distinction to the "fat" random 
graphs that appear in the planar Feynman diagram limit, N — » oo, more familiar from discretized 
models of two dimensional gravity. 

The interest of the thin graphs is that they give mean field theory behaviour for spin models living 
on them without infinite range interactions or the boundary problems of genuine tree-like structures 
such as the Bethe lattice, g-state Potts models display a first order transition in the mean field 
for q > 2, so the thin graph Potts models provide a useful test case for exploring discontinuous 
transitions in mean field theories in which many quantities can be calculated explicitly in the saddle 
point approximation. 

Such discontinuous transitions also appear in multiple Ising models on thin graphs and may have 
implications for the use of the replica trick in spin glass models on random graphs. 



1 Introduction and Reprise of Continuous Transitions 



A simple and elegant method of describing spin models on random graphs, drawing inspiration from 
the matrix model methods [Q used to describe planar random graphs in two-dimensional gravity, was 
first proposed in |2). It was observed that the requisite ensemble of random graphs of unrestricted 
topology could be thought of as arising from the perturbative Feynman diagram expansion of a scalar 
integral in much the same manner as the planar graphs that appear in two dimensional gravity theories 
were generated from the perturbative expansion of a matrix integral. In effect, the unrestricted random 
graphs appear in the N — ► 1 limit of an N x N Hermitian matrix model, which we denote as "thin" 
graphs, to distinguish them from the planar "fat" graphs which appear in the N — > oo limit and still 
retain their matrix structure. Throughout the paper we will use "thin graphs" and "Feynman diagrams" 
interchangeably to denote the random graphs of unrestricted topology on which our spin models live. 
Spin models on random graphs are of interest as they will display mean field behaviour because the 
graphs have a tree-like local structure The advantage of using random graphs, which are closed, 
over genuine tree-like structures such as the Bethe lattice is that dominant boundary effects are absent. 
The complications, both analytical and numerical, of being forced to consider only sites deep within the 
lattice are thus absent. Other ways of accessing mean field behaviour, such as infinite range interactions, 
are not very well suited for numerical simulation. 

In previous papers we showed that the thin graphs of the Feynman diagram expansion offered a 
practical method of investigating mean field models both analytically and numerically The equilibrium 
behaviour of ferromagnetic Ising models j| and spin glasses H was found to parallel that of the equivalent 
model on the appropriate Bethe lattice with the same number of neighbours and the analytical 
treatment offered a different perspective to previous approaches to random graph spin models and spin 
glasses 0, ||] . The investigation of dynamical phenomena such as aging effects in spin glasses || was also 
found to be facilitated by random graph simulations. Other authors have also employed random graphs 
in simulations of the random field Ising model |To[ in order to avoid boundary problems with the Bethe 
lattice. 

Analytical calculations using the approach of Q involve simple saddle point methods for standard 
integrals, or quantum mechanical path integrals in the case of continuous spins [ jlT| . If we consider 
undecorated random graphs, taking a 3 theory for definiteness which will generate 3-regular random 
graphs [], the number of such graphs with 2n vertices can be calculated as 

which, when evaluated using a saddle point approximation, gives the correct counting 

m 2 " (6n-l)!! 
" \6J (2n)H ' [) 

To include an Ising model we now decorate the vertices of the graphs with Ising spins having a Hamiltonian 

<ij> 

where the sum is over nearest neighbour sites. The partition function is then given by 
where K is defined by 

*2 = ( \ 7 ) < 5 > 



1 We will restrict ourselves to 3 or 3-regular random graphs throughout. The saddle point equations may still be solved 
with larger numbers of neighbours, but become rapidly more complicated. 



1 



and the action itself is a direct transcription of the matrix model action [T^| to simple scalar variables 

S = \ E ^ K ab<Pb ~ 5 (4 + <t ) ■ (6) 

a, b 

The sum in the above runs over ± indices The coupling c = exp(— 2/3) and the + field can be thought 
of as representing "up" spins with the 0_ field representing "down" spins. It is necessary to include the 
counting factor N n to disentangle the factorial growth of the undecorated graphs from any non-analyticity 
due to phase transitions in the decorating spins. One is also obliged to pick out the 2n-th order in the 
expansion explicitly with the contour integral over A as, unlike the planar graphs of two dimensional 
gravity, A cannot be tuned to a critical value to cause a divergence. 

The mean field Ising transition manifests itself in this formalism as an exchange of dominant saddle 
points. Solving the saddle point equations dS/d<f)± = 

4>+ = 4>\ + c</>- 

0- = 2 _+c0+ (7) 

(which we have rescaled to remove A and an irrelevant overall factor) we find a symmetric high temperature 
solution 

4>+ = <t>- = 1 - c (8) 
which bifurcates at c = 1/3 to the low temperature solutions 



1 + c + \/l - 2c - 3c 2 



1 + c - VI - 2c - 3. 
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(9) 

The bifurcation point is determined by the value of c at which the high and low temperature solutions 
for are identical, which appears at the zero of the Hessian det(d 2 S/d(f> 2 ). The magnetisation order 



parameter for the Ising model can also be transcribed directly from the matrix model 1 12 



4>l -0 3 

M = ^± — ^ 10 

03_ + 0i v ' 

and shows a continuous transition with mean field critical exponent {(3 — 1/2). The other critical 
exponents may also be calculated and take on mean field values. 

Simulations of the Ising model are in very good agreement with the analytical results Q even on 
a single graph, which at first sight is rather surprising as the saddle point calculations are formally 
for an annealed ensemble of graphs. This appears to be true for all models where one might expect 
self-averaging, such as ferromagnetic Ising and Potts models - one large graph can be thought of as a 
collection of smaller graphs in these cases. With spin glasses it is still obligatory to consider a quenched 
ensemble of random graphs in order to take a (quenched) average over the disorder. 



2 Potts Models 

The Hamiltonian for a q-state Potts model can be written 

h = pJ2( s °^- 1 ) ( u ) 



<y> 



where the spins Oi take on q values. The matrix model actions for such Potts models are well known 
(though only solved exactly so far for q = 3 For the 3-state Potts model the action is 



>' = l&l + 4>l + 4>l) - c(0i0 2 + 0i03 + 2 3 ) - + 4>l + 01), ( :i->) 



We have rescaled the 4>'s with respect to W, Bl tj for uniformity with the Potts model notation. 
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which can be used as the action on thin graphs if one takes, as in the Ising case, the 0's to be scalar 
variables. For a q-state Potts model c = l/(exp(2/3) + q — 2). 

"Ising- like" solutions to the 3 and 4-state Potts models were presented in M , the 3 state case being 

01,2,3 = l-2c (HighT) 

(LowT), 



1 + Vl - 4c - 4c 2 

'1,2 = n 



1 + 2c - Vl - 4c - 4c 2 



(13) 



These high temperature and low temperature solutions are equal at the zero of the Hessian c = 1/5, i.e. 
g = 4. Similarly, the 4-state Potts model has the action 

S = \{4>\ + 02 + 03 + <t>l) ~ C(0102 + 0103 + 0104 + 0203 + 0204 + 030 4 ) " ^ (01 + 02 + 03 + $)> ( 14 ) 

and solving the saddle point equations again gave Ising-like solutions 

01,2,3,4 = 1 — 3c {HighT) 



'1,2,3 



1 - c + \/l - 6c- 3c 2 



1 + 3c - Vl - 6c - 3c 2 



(LowT), 



04 = — 2 (15) 

where the solutions matched at c = 1/7, g = 5 []. The picture is repeated for higher q, where the action 
is 

z—1 i<j i—1 

and one finds a high temperature solution of the form 0j = 1 — (g — l)c,Vi bifurcating to a broken 
symmetry solution 0j = . . . g _i 7^ 9 at g = q + 1. 

These results are somewhat surprising on two counts. Firstly, the motivation for using thin graphs 
was that they provided easy access to mean field results. However, it is known that the mean field theory 
for Potts models predicts a first order transition for q > 2. Secondly, all the thin graph results so far 
for various models have been identical to the corresponding Bcthc lattices, even down to non-universal 
features like the transition temperatures. Explicit calculations on the Bethe lattice have also given first 



order behaviour for q > 2 14 and shown the values of g obtained above for the 3 and 4 state models 
correspond to spinodal points on the Bethe lattice. The models hit a first order transition before attaining 
these points. One might therefore expect that a first order transition should be lurking in the saddle 
point solutions for the actions above, given the previous tendency for the thin graph results to slavishly 
parallel the Bethe lattice. 

The resolution of the conundrum is implicit already in the solutions in equs.([l3],[l5]) and their higher 
q equivalents. If we look at the Ising solution of equ.(^) we can see that the low temperature branches 
become real exactly at the transition point, whereas the square roots in the q > 2 Potts solutions become 
real at larger c, and hence higher temperature. The topology of the phase diagram is perhaps best 
understood by plotting the magnetisation against c, which we do in Fig.l for the Ising solution and in 
Fig. 2 for the 4 state Potts model (all other q > 2 state model s being of similar form). Low temperature 
corresponds to small c and high temperatures to large c for all the Ising and Potts models. For conformity 
with the Potts notation we define a Potts style magnetisation for the Ising model as 

= + 0- ri7 x 

(0^ + 03_)' (0^ + 03) {l) 



3 The right hand side of the LowT solutions may be exchanged, as oije might have expected on symmetry grounds, and 
we have used this freedom to put the solutions in a tidier form than in 
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where the first variant gives the upper (solid) branch of the low temperature magnetisation curve in Fig.l 
and the second the lower (dashed) branch. Both definitions give the same value for the horizontal (dotted) 
high temperature solution with m — 1/2 as <f) + = 4>- there. We can see that the paramagnetic high 
temperature solution bifurcates at O (c=l/3) to give the upper and lower magnetised low temperature 
branches. 

The equivalent magnetisation m for the Potts models is defined as 



on the upper low temperature branch where it gives a maximum value and 



(18) 



(19) 



on the lower low temperature branch, where it gives a minimum. On the symmetric high temperature 
branch (dotted) both the definitions are, of course, equivalent and one finds m = X/q. The standard 
Potts model order parameter can then be defined as 

M=(? max(m)-1 
1 - 1 

which is zero in the high temperature paramagnetic phase and tends to one in the magnetised low 
temperature phase. 

The solutions of equs.(|l5|) give the lower (dashed) branch in Fig. 2, descending from the point O at 
which the square roots become real. We can see that the magnetisation "pitchfork" of the Ising diagram 
becomes skewed for q = 4 (and all other q > 2). The most important feature is that the upper (solid) 
branch does not connect continuously with the (dotted) high temperature solution, which joins the lower 
branch at P. The upper branch is simply obtained for all q by choosing the opposite sign for the square 
roots to the lower branch. As can be seen in Fig. 2 it corresponds to the true low temperature magnetised 
phase where m — > 1 (and hence M — > 1) as T — ► 0. The first order transition, denoted by a vertical line 
with ends labelled Q in Fig. 2 takes place when the free energy of the upper branch is equal to the free 
energy of the high temperature solution. In the saddle point approximation the free energy to lowest 
order in the number of vertices n is just the action S so the first order transition point is given by the c 
value, and hence temperature satisfying 

S{upper branch) — S(high temperature). (21) 

As one can see in Fig. 2 as the temperature (i.e. c) is reduced a first order transition intervenes between 
the Q's before P is reached. Similarly, as the temperature is increased from zero along the upper branch 
a first order transition occurs before O is reached. The portions of the magnetisation curve PO and the 
dotted horizontal line to the left of P represent unstable states, whereas PQ, QO and the lower dashed 
branch to the left of P represent meta-stable states The meta-stable portions of the curve would be 
accessible by superheating out of the magnetised phase (QO) or supercooling from the paramagnetic 
phase (QP, P — -> origin). For completeness, we have listed the stable low temperature solutions for 
q = 3, 4, 5, 6 state Potts models in Appendix. A, which are the q values simulated in the next section. 

Various features of the solutions merit discussion. Equ.(|2l|) can be solved analytically for moderate q 
values on </> 3 graphs without too much difficulty, and all the solutions fit the following compact formula 
for the critical value of c at Q 

1 - fa - I)- 1 / 3 

c(Q) = lg _o ■ ( 22 ) 

q z 

Indeed, if one takes the conjectured q state solutions in Appendix. A it is possible to write down the saddle 
point action on the (upper) low temperature branch in terms of <f> = 4>\... q -i and = (j) q by substituting 

4 These considerations, and indeed the two figures, are essentially identical to those for the Bethe lattice in JmJ. 
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them into equ. (|l6|) 



S= l -( q -l)[l-c{ q -2)]cj ) 2 - 

and similarly on the high temperature branch 

Q 



5<" 



1 TO 



1 ~ 



s 



(1- C (g_l))^-10; 



(23) 



(24) 



where 



= 1 — (q — l)c. Setting (equ.(|23|)) and (equ.(p4|)) equal, one also obtains the same c{Q) as in 
equ. (22) above. 

It is possible to calculate the jump in the magnetisation AM along the vertical line at Q and one 
finds in all cases 

AM = (25) 

q - 1 

If one now refers back to the Bethe lattice calculations of jHj one can see that AM is identical to that 
observed on the Bethe lattice. In addition, allowing for the differences in conventions ^, the formulae for 
the critical coupling c(P) in equ. (22) is also identical to that for the Bethe lattice transition. The zero of 
the Hessian for the q state Potts model action gives us the value of c at P where the high temperature 
solution joins the lower branch 

<P) = ^ (26) 

and if we assume that the conjectured q state low temperature solution in Appendix. A is correct, we can 
also calculate the value of c at O where the square roots become real 



c(0) 



q- 1 -2y^T 
(?-!)(? -5) 



(27) 



(q = 5 can be handled by taking the limit q — > 5). As q is increased the separation between points points 
O and P increases. In Table. 1 below we list for convenience the c values of points O, P and Q for the 
(7 = 3,4,5,6 state Potts models. 



q 


c{0) 


c(P) 


c(Q) 


3 


0.20711 


0.20000 


0.20630 


4 


0.15470 


0.14286 


0.15332 


5 


0.125 


0.11111 


0.12335 


6 


0.10557 


0.09091 


0.10380 



Table 1: The c values for the points O, P, Q along with j3 C rit for q — 3, 4, 5, 6 state Potts models 
The values of c{0) and c(P) that we have found are again identical to those on the Bethe lattice. 



3 Simulations 

The acid test of the saddle point solutions is whether they match up with simulations. We do not 
attempt a high accuracy verification of the analytical results here, but rather a consistency check on the 
first order nature of the transition and a verification of the values for c(Q) and AM calculated in the 
previous section. To this end we generated single </> 3 graphs with 250, 1000 and 2500 vertices for each of 
the q state models. We verified the results by repeating all the simulations on a second, different graph 
for each size with identical results within the error bars in all cases. 

5 The 9 of |14| is equal to exp(— 2/3) in our notation. 
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The generation of the graphs is easier than the planar <fr 3 graphs used in 2c? gravity simulations 
because of the absence of a constraint on the topology. This obviates the need to use (for instance) 
Tutte's algorithm in producing the graphs. The simulation itself used the Wolff algorithm, which will 
not beat the super-exponential slowing down right at the first order transition point, but is expected to 
be efficient elsewhere. We simulated a range of (3 values, allowing 20,000 equilibration sweeps followed 
by 20,000 xJV cluster updates, where N was O(10) - the exact value depending on the mean cluster size. 
Measurements were made every N cluster updates of all the standard thermodynamic quantities, the 
energy E, the magnetisation M, specific heat C and magnetic susceptibility X- We also measured various 
Binder's cumulants for the magnetisation < M 4 > / < M 2 > 2 and energy < E 4 > / < E 2 > 2 as well as 
correlation functions and autocorrelations. 

We focus our attention first on the magnetisation curves for the various models, which are presented, 
with fortuitous numbering, in Figs. 3, 4 for the q = 3,4 state models on graphs with 250,1000 and 2500 
vertices. The smaller graphs display greater finite size rounding, but by the time one has got to 2500 
vertices the agreement with the magnetisation calculated from the saddle point solutions (which are 
formally for an infinite number of vertices) is already quite good . On both plots we have delineated 
the expected critical points and jumps in the magnetisations. This agreement deteriorates somewhat 
for a given lattice size as q and hence the strength of the transition increases. The j3 crit for various 
q are estimated from the simulations by looking at the crossing of Binder's magnetisation cumulant 
< M 4 > / < M 2 > 2 for the various graph sizes. As one can see in Fig. 5 for the 3 state Potts model 
(which is representative) the errors in the measurement of the cumulant are quite large, but even given this 
the estimated critical temperatures are all close to those calculated in the saddle point approximation. 
We also list the (3 values for the two spinodal points O and P for comparison. As one can see even 
the rather modest simulations carried out here are sufficient to show that the spinodal point P can be 
excluded as the transition point in all cases. The results for all but the 3 state model cannot definitively 
exclude the other spinodal point O as the critical point, but the first order nature of the transition and 
the value of the jump in the magnetisation, as discussed below, favour a transition at Q as predicted by 
the saddle point calculations. 



q 


3 


4 


5 


6 


ftcrit 


0.674(2) 


0.75(1) 


0.81(1) 


0.87(1) 


P{Q) 


0.67369 


0.75451 


0.81533 


0.86441 


(3(0) 


0.67122 


0.74804 


0.80472 


0.84986 


m 


0.69315 


0.80470 


0.89588 


0.97295 



Table 2: The estimated /3 C rit, along with 
the calculated (3(Q) , (3{0) , (3{P) . 

The AM values are estimated by eye-balling the magnetisation curves for largest graphs and are 
consequently to be taken with a larger pinch of salt than the other measurements, but they are all 
consistent with the (q — 2)/(q — 1) calculated in the previous section. We tabulate the results for the 
magnetisation jumps measured from the simulations below in Table. 3 along with (q — 2)/(q — 1) for 
comparison. 



q 


3 


4 


5 


6 


AM 


0.50(5) 


0.68(1) 


0.72(5) 


0.9(2) 


(«-2)/(g-l) 


0.5 


0.66667 


0.75 


0.8 



Table 3: The estimated AM along with 
the calculated (q — 2)/{q — 1). 

Further confirmation that the transitions are indeed first order can be obtained by looking at the 
values of the Binder's energy cumulant near (3 cr it- This is expected to scale to 2/3 for a continuous 
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transition and a value less than 2/3 if the transition is first order. The values are listed in Table. 4 and 
show a clear tendency to decrease around fl C rtt that grows stronger with increasing q. 



q 


3 


4 


5 


6 


<E 4 > 
<E 2 > 2 


0.664 


0.65 


0.64 


0.60 



Table 4: The values of Binder's energy cumulant at the estimated f3 cr u 

The energy itself is calculable from E = —d\ogZ/d(3 and is discontinuous at a first order transition. 
One finds values that are again in good agreement with the simulations by the time one has reached 2500 
vertices. The measurements for the 4 state Potts model are shown in Fig. 6. In a similar vein quantities 
such as the specific heat C = 2 {d 2 \ogZ/d0 2 ) or the magnetic susceptibility may be calculated from 
the saddle point solutions and all give very satisfactory agreement with the measured quantities in the 
simulations. We do not describe these here as it is clear from the results already presented that, even 
given the limitations of the fairly modest simulations, there is ample support for the correctness of the 
saddle point solutions and the picture of the first order transition that they suggest. 

In closing, we note that something akin to a standard finite size scaling analysis is possible with the 
thin graph approach to simulations, as witnessed by the use of the Binder's cumulant to estimate the 
critical temperature in the current work. The place of the factor L~^' v that appears in finite size scaling 
on standard lattices, where L is the linear size of the lattice, is taken by rT x l vd where n is the number 
of vertices in the graph. Although d is formally infinite, the combination vd is still well defined and all 
the scaling relations may be written in terms of this. Entirely analogous tactics have been used in the 
analysis of simulations of spin models in planar diagrams in theories of two-dimensional gravity, where d 
in this case was a dynamically generated fractal dimension that was a priori unknown. On the analytical 
side 1/n corrections may be obtained (and have been obtained already for the Ising model in Q) by 
calculating the determinantal corrections to the saddle point solutions, so corrections to scaling can be 
obtained. Such issues would be worth pursuing if very high accuracy verification of the correspondence 
between calculation and simulation were required. 



4 Discussion 

Our previous analytical and numerical work Q [s[ ^| on spin models on Feynman diagrams had concen- 
trated on the case of continuous transitions. The results in the current paper show that the simple saddle 
point equations that determine the phase structure of such models are also adequate to describe first 
order transitions. In the continuous case the critical point was pinpointed by finding the zeroes of the 
Hessian and corresponded to a bifurcation of magnetised states from the unmagnetised high temperature 
solution. This gives one the lower temperature spinodal point P in the Potts models, the true first order 
transition point at Q being determined by matching the saddle point actions (i.e. free energies) on the 
two branches of the solution. The upper spinodal point O is fixed as the point at which a square root 
appearing in the low temperature branches becomes real. 

The critical temperatures we have calculated, the jump in the magnetisation and the magnetisation 
curves themselves are identical to the results obtained in |L4j on the Bethe lattice. We thus conclude that 
the (mostly large) loops that are present in the Feynman graphs have no effect of the critical behaviour 
of ferromagnetic Potts models by comparison with the corresponding Bethe lattices. This is consistent 
with the earlier work on continuous transition which also demonstrated Bethe-lattice-like results. The 
loops will, however, have an effect in the antiferromagnetic models considered in EJ , where a two-step 
invariant measure which presupposes a bipartite lattice was instrumental in the solution. As already 
noted in |I| for the antiferromagnetic Ising model, loops of both even and odd length are present in the 
Feynman diagrams, so frustration will be present on 4> z graphs if q < z. This apparently leads to a spin 
glass phase rather than antiferromagnetic order. In matrix models it is possible to arrange only even 
sided polygonations by using complex rather than Hermitian matrices, but it is not clear to us how to 
perform a similar trick on generic - rather than planar - Feynman diagrams. 
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We find it rather remarkable that it is possible to write down the actions in equs.( |2^f24| ) for both 
the high and low temperature branches for arbitrary q, based on the ansatz for the q-state solutions in 
the Appendix. This is on a par with the results in || for the Hessian of an arbitrary number of Ising 
replicas on thin graphs. The availability of a solution where q appears explicitly as a parameter opens 
the possibility of exploring the q — > 1 limit of the model, related to percolation, which we will address 
in a further publication. The percolative transition on the Bethe lattice has some unusual features [jDsl , 
which one might also expect to manifest themselves on thin graphs. The behaviour of the model in an 
external field, considered in some detail on the Bethe lattice in |Q, can also be investigated with the 
thin graph formalism both analytically and numerically without much by way of complications over the 
results presented here. 

First order transitions also appear in the so-called Ising replica magnet , the finite k version of the 
fc-Ising replicas used in taking the k — ► limit ("replica trick") with quenched disorder. In these models 
on thin graphs the spin glass transition temperature appears as a sort of spinodal point || for all k > 2, 
with the transition being continuous for k — 2. A first order (in the overlap) transition intervenes before 
the putative continuous spin glass transition is reached for k > 2. If it were not for this one would be 
tempted to argue that the k — > limit for the spin glass transition temperature was trivial, as it is the 
same for all k. The role of the first order transition to a replica symmetric state when k > 2 and the 
nature of the k — > limit for thin graph models thus requires further elucidation. It is possible that the 
Potts models results presented here may cast some light on its properties. 

The numerical work in this paper was intended as a consistency check of the formalism, rather than 
a full scale numerical investigation and finite size scaling analysis of the models. Nonetheless, it is clear 
from the results presented that the transitions are first order as predicted. The agreement between the 
simulations and the saddle point calculations for the critical temperatures and the observed jumps in the 
magnetisations are very satisfactory even on graphs with 2500 vertices. 

In summary, we have seen that the thin graph approach is a convenient way of performing calculations 
and simulations for mean field Potts models with first order phase transitions. Such models, as well as 
being of interest in their own right, may help in understanding mean field spin glass and percolative 
transitions. 
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Appendix A 

We list here the low temperature solutions (upper branch in Fig. 2) for the q = 3, 4, 5, 6 state Potts models 
which arc referred to in the text. Note that the solutions appear follow a regular pattern and we have 
indicated the conjectured q state solution (which fits the cases listed below, the Ising model q = 2, and 
all the other higher q solutions we checked explicitly before exhaustion set in) at the end. Just as for the 
lower branch solutions where the signs in front of the square roots are reversed, the right hand sides of 
the solutions can be exchanged for a given q. 

3-State 



1 - VI - 4c - 4c 2 
01,2 = 2 

1 + 2c + Vl - 4c - 4c 2 



4- State 

1 - c - VI - 6c - 3c 2 

01,2,3 = ^ 

l + 3c+Vl-6c-3c 2 
04 - g 

5- State 

1 - 2c - VT - ^ 
2 

1 + 4c + VI - 8c 



'1,2,3,4 



»b = 



6-State 

1 - 3c - V5c 2 - 10c + 1 

01,2,3,4,5 — ~ 



1 + 5c + V5c 2 - 10c + 1 

»6 = » 



q-State 

1 - (g - 3)c - y/1 - 2(g - l)c + (g - 5)(g - l)c 2 

01.. .9-1 = o 



1 + (g - l)c + v/l-2(g-l)c+(g-5)(g-l)c 2 
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Figure 1: The magnetisation m for the Ising model as calculated from the saddle point solutions. The 
high temperature branch is shown dotted, the upper low temperature branch solid and the lower low 
temperature branch dashed. 



11 



1 




Figure 2: The magnetisation m for a 4 state Potts model as calculated from the saddle point solutions. 
The linestyles are as in Fig.l and only the portion of the graph close to the transition point is shown for 
clarity. 



12 



i r 



M, 3 state Potts 



0.5 



X^ X 

M X 



x 



X 



<> 



X 



X 



X 



o o I I 



X 



o 



o x 




0.6 



JJ L 



0.65 



0.7 
(3 



x 2500 
o 1000 
x 250 



0.75 



0.8 



Figure 3: The magnetisation M (— (qmax(m) — l)/(q — 1)) for a 3 state Potts model as measured in 
the simulations. The saddle point solution is shown as a dashed line. The horizontal bar represents the 
height of the "jump" — ► AM in the magnetisation. The error bars are too small to be seen in all but 
the central points 
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Figure 4: The magnetisation M for a 4 state Potts model as measured in the simulations. Key as for 
Fig.3. 
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Figure 5: The crossing of Binder's magnetisation cumulant for the 3 state Potts model as measured in 
the simulations. 
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Figure 6: The energy for a 4 state Potts model as measured in the simulations. The position of the 
transition as calculated in the saddle point approximation is shown again as a dotted line. 
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